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ABSTRACT 

The  application  of  a  simple  discrete  velocity  model  to  lov  Mach  number 
Couette  and  Rayleigh  flow  is  investigated.  In  the  model,  the  molecular  ve¬ 
locities  are  restricted  to  a  finite  set  and  in  this  study  only  eight  equal 
speed  velocities  are  allowed.  The  Boltzmann  equation  is  reduced  by  this 
approximation  to  a  set  of  coupled  differential  equations  which  are  shown  to 
be  identical  in  form  to  those  produced  when  the  seme  approximation  is  applied 
to  the  Krook  equation.  The  fluid  velocity  and  shear  stress  in  Couette  flow 
are  in  approximate  accord  with  those  of  Vang  Chang  and  Uhlenbeck  and  of  Lees 
over  the  complete  range  of  Khudsen  number. 


( 


INTRODUCTION 

Is  a  previous  report  (l),  an  approximate  description  of  lov  Mach  number 
rarefied  shear  flow  vas  obtained  by  restricting  the  molecular  velocities  to 
a  finite,  prescribed  set,  an  approximation  which  reduces  the  governing  lntegro- 
differential  equation  to  a  set  of  coupled  differential  equations.  That  study 
vas  based  on  the  Krook  approximation  to  the  Boltzmann  equation  and  eight,  equal 
magnitude,  molecular  velocities  were  allowed.  It  is  the  purpose  of  this 
report  to  show  that  when  the  same  approximation  is  applied  to  the  Boltzmann 
equation  itaelf,  the  resulting  set  of  differential  equations  is  identical 
in  form  to  that  obtained  from  the  Krook  model;  the  only  difference  is  in  the 
magnitude  of  the  relaxation  time.  In  the  Krook  model  this  parameter  vas  set 
equal  to  the  mean  time  between  collisions,  x/f;  here  it  is  shown  that  it  is 
approximately  X/.68  c.  The  Couette  flow  results  of  Ref.  1  are  brought  into 
approximate  quantitative  agreement  with  those  of  Wang  Chang  and  Uhlenbeck  (2) 
and  of  Lees  (3)  by  the  use  of  the  new  relaxation  time.  There  is  also  significant 
Improvement  in  the  Rayleigh  flow  solutions. 

THE  DISCRETE  VELOCITY  MODEL 

The  Boltzmann  equation  may  be  written 


df 

TSt 


a) 


in  which  f  is  the  distribution  function  depending  on  the  time  and  space 
variables  t,  x,  y,  i  and  on  the  molecular  velocity with  components  u,  v,  v. 
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Tt)9  rate  of  change  of  f  due  to  collisions  Is  written  as  tbs  gain  sinus  the 
loss,  0  -  L.  Bov  if  the  molecules  occupy,  between  collisions,  a  finite  set 
of  cells  located  at  ^  in  velooity  space,  (l)  can  be  written 


in  which  Hi  is  the  niaber  of  nolecules  per  unit  voluse  with  velocity  v±. 

Before  evaluating  (Qj-I^)  for  eight  cells,  consider  as  a  simpler  example 
the  two  -  dlaensional  gas  lying  in  the  u,  v  plane  with  four  velocities  as 
sketched  in  Fig.  la. 


Figure  1 
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The  speed  la  taken  to  be  the  mean  thermal  speed,  e.  The  rate  of  change  of 
Nj,  for  Instance,  due  to  collisions  is  determined  as  follows.  Loss  of  Mole¬ 
cules  from  cell  1  occurs  only  when  these  molecules  collide  with  those  in 
cell  4,  for  collision  with  cell  2  or  3  occupants  results  only  in  ah  exchange 
of  cells  (see  Figs,  lb  and  lc).  Thus  the  loss  from  cell  1  can  be  written 

where  v  is  the  relative  speed,  equal  to  2  c,  and  8  is  the  effective  col- 

r  8 

llslon  cross  section,  l.e.,  the  cross  section  which  deflects  collision 
partners  from  cells  1  and  4  to  2  and  3.  These  same  arguments  show,  of 
course,  that  molecules  are  thrown  into  cell  1  from  collisions  2  -  3i  there¬ 
fore 

0X  -  2  c  Se  H2  »3 

and 

(sr)c-2'8.  (»2 >3 -  «! \>- 

The  expressions  for  the  other  cells  can  be  written  in  the  same  way* 

When  there  are  eight  cells,  synsetrlcally  placed  in  the  eight  velocity 
quadrants  (see  Fig.  2),  and  the  flow  Is  Independent  of  z  and  symmetrical 
about  the  w  «  0  plane, 

-  »5,  »2  -  H6,  k3  -  iy  ^  -  H8 


> 


V 


Figure  2 

Under  these  conditions  and  with  the  assumption  that  the  molecules  are  hard 
elastic  spheres,  it  Is  shown  In  the  Appendix  that 

(sr)  “  (1  +  ^  ®  8  <*2  N3  -  *1  V  (3) 

c 

where  8  Is  the  mutual  collision  cross  section. 

Anticipating  the  application  of  these  equations  to  low  Mush  nvnber 
Couette  and  Rayleigh  flow,  let  us  divide  Iq.  (2)  by  the  masher  density,  N, 
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assuned  constant .  Than  making  use  of  Eq.  (3)  and  corresponding  expressions 
for  the  other  cells  ve  get. 


dn.,  in.  dn^ 

W  '  4  ST  +  4  ST 


2  9  (n2  n3 


(U) 


8n2  dn2  an« 

sr  +  4*r +  4sr 


2  9  (nx  n^ 


■W 


(kb) 


Sn,  Sno  .  x 

q  ^-A  -  q  ■  2  9  (n^  n^  -  ng  n^) 


3ni  3n.  dni.  .  ... 

sr*4sr"  4sr-2e  <°2”3‘nink)  144 

where  n^  is  the  fraction  of  molecules  in  cell  i,  q  *  c/JT>  and 

0  -  pc  (  1  +  2//3)  c  S  N. 

v  6 

Several  general  features  of  these  equations  are  worth  noting.  First, 
when  they  are  multiplied  successively  hy  the  colliaional  invariants,  m, 
mu  j,BVj  (ils  the  molecular  weight)  and  sunned,  the  right  hand  sides 
vanish,  a  consequence  of  the  fact  that  the  collisions  sketched  in  Fig.  1 
conserve  the  number  of  molecules  and  satisfy  the  equations  of  motion. 
(Energy  is  autcnaticsdly  conserved  since  the  speed  is  constant.) 
lfce  equilibrium  condition  is 

ni  n4  “  n2  n3 
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a  remnant  of  tba  general  equilibrium  condition 

t(7'  t&J  Kfj,). 

If  the  n^  'a  are  Independent  of  x  and  y  it  la  eaay  to  ahov  with  the  help  of 
Bqa.  (4)  that  the  H  function,  defined  by 

H  "  I  ni  ^  V 
1 

obey*  the  equation 

dt  "  2  *  (a2  a3  "  al  V  (la  al  nk  “  1X1  a2  n 3 ) 

and  hence 

dt 

aa  In  the  exact  equation 

COUTTTK  AKD  RATUtlOH  TUN  EQUATIONS 

To  apply  Eqa.  (k)  to  Couette  and  Rayleigh  flow  we  aaatane,  a a  waa  done 

In  Ref.  1,  that  the  fluid  velocity  In  the  y  direction,  V,  la  negligible 

(the  walla  are  parallel  to  the  x  -  aids)  when  the  lfcch  number  la  low  enough. 

Then  aince  V  »  I  n,  v., 

1  1  1 

nl+D5+n2+n6"n3+n7+nfc+n8*2 
and  fran  the  symsetry  about  w  -  0, 

”i  *  *  "3  *  "k  ‘  I- 
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These  relations  make  Rqs.  (4b)  and  4d),  say,  superfluous  and  allow  the  others 
to  be  written 

an  an. 

at1  +  4  "Sy  “  2  (  -  nl  +  V  (5a) 

an,  an,  a 

5t  •  4  3y  “  ■  2  (  ‘  nl  +  n3)  (5b) 

when  the  terns  depending  on  x  are  dropped.  These  are  Bqs.  (l8a)  and  (l8c) 

of  Ref.  1  with  0  appearing  in  place  of  a ,  quantities  which  differ  only  by 
* 

a  numerical  factor. 

This  Is  the  result  stated  above,  namely,  that  to  the  present  approxi¬ 
mation  the  Boltzmann  and  KTook  equations  lead  to  differential  equations 
describing  Couette  and  Rayleigh  flow  which  are  identical  In  form.  This 
could  perhaps  be  taken  as  simply  an  Indication  of  the  crudity  of  the  eight 
cell  model  were  It  not  for  the  many  points  of  agreement,  discussed  in  the 
following,  between  the  results  given  by  this  simple  model  and  those  of  much 
more  rigorous  analyses. 

In  the  Krook  model  the  magnitude  of  the  effective  collision  rate  must 
be  chosen  arbitrarily  -  It  Is  of  the  order  of  the  Inverse  of  the  mean  time 
between  collisions  -  and  In  Ref.  1  It  was  set  equal  to  this  quantity,  l.e., 

*  In  Ref.  1,  n,  represents  the  sum  of  the  fraction  of  molecules  In  cells 
1  and  5#  l.e.,  2  n^  represents  2  ng;  etc.  Iqp.  (5),  being  linear,  can 
be  similarly  Interpreted. 


to  o  -  c/X.  It  li  shown  in  ths  Appendix  thet  8  -  ,68  c/X.  With  this  de¬ 
gree  of  arbitrariness  removed  from  the  discrete  velocity  aodel,  it  is 
worthwhile  to  ccnpere  the  solutions,  quantitatively,  with  those  of  Hang 
Chang  and  Uhlenbeck  (2)  and  of  Lees  (3).  Only  Couette  flow  will  be  dis¬ 
cussed  here;  the  Reyleigh  flow  solutions,  described  in  Ref.  1.,  are  also 
inproved  when  8  replaces  o  but  the  details  will  be  given  in  a  later  report. 
COUKTTX  FLOW  SOLUTION 


Figure  3. 


Consider  the  steady  Couette  flow  sketched  in  Fig.  3.  For  this  problem 
Equations  (5a)  and  (5b)  can  be  written 
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d  n_ 

- s  -  Qf(  -  n.  +  n.) 

d  y  3 


(6a) 


-1  -  «(-  n.  +  n,)  (6b) 

d  y  3 

vhere  y*  ■  y/d  and  a  «  d  6/2  q 

When  the  reflection  from  both  walla  is  diffuse  and  n^  -  n^,  Ug  -  ng, 
etc.,  the  boundary  conditions  are 


vhere  Uv  is  the  difference  in  the  wall  velocities.  These  equations  state 
that  the  average  x  -  component  of  velocity  of  the  molecules  leaving  the 
upper  vail  is  -  U J2  at  the  upper  wall  and  that  of  the  upward  flowing 
particles  is  U J2  at  the  lower  wall.  Then  since  n^  +  n2  ■  n^  +  n^  ■  j-, 

*>]_(-  “  $0-  "  V2q)*  n3^  *  +  V21) 

Solutions  to  Equations  (6a)  and  (6b)  satisfying  these  conditions  sre 
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y*  + 1  - 


2(or  +  1}  q 


■i-i[T5TTr» 

5[vTi)  f  y*  *  1  *  Z&Vfl  i] 


°3'3 

Frco  these  'distribution  functions '  the  fluid  velocity  is  the  r 
direction,  U,  le  found  to  he 


U  -  £  ut  nA  -  2  q  [  -  nx  +  n2 

Of  * 

"  •  T«hTT  uv  y 

and  the  shear  stress,  p  ,  is  given  hy 


“s  *  M 


»  [i  -  M»!  *  >3)] 
(7) 


"  l  \  \  ”1  •  2  "  «*  [  ■  °1  *  ”2  *  ”3  •  *]  ■  11  P  *Vl  *  °3 

<8) 

Defining  the  viscosity,  p,  hy  p^/  ve  have  frcn  Equations  (7)  end  (8)$ 

U  -  P  q  d/2  o  -  p  q2/» 
and  with  0  -  .68  c/X  and  q  ■  c/fe 
U-  .49  P  c  X 

coopered  with  the  classical  value  of  .499  P  c  X. 

Sext  the  fractional  slip  velocity,  2  A  o/u^,  defined  hy 
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2  A  U/Uy 


Uv/2'U 

^T5 


til 


is,  from  Equation  (7), 

2  2  "/"»  ■  s4t-  ■  *  1  <9) 

Lees  finds,  for  Maxwell  molecules, 

2  2  •  -T^m 

These  expressions  are  shown  in  Figure  k ,  together  with  the  results  of  Mang 
Chang  and  Ublenbeck  who  also  treated  Maxwell  molecules . 

Finally,  the  ratio  of  the  shear  stress  to  the  free  molecule  value, ? 
pyx/(pyx)f  n  i  can  he  found  from  Equation  (8)  to  he  given,  also,  hy 

vyx/(vyx*t.a.  "  cr  +  1  “  .59  d/X  +  1 

Similarly,  Lees'  expression  is 

pyx^pyx^fds.  “  .5  d/x  +  1 

In  Figure  k  the  Hang  Chang  and  Uhlenbeck  values  are  compared  with  these  two 
equations.  Equation  (10)  lies  helow  those  of  the  other  investigators  at 
large  d/x  because  the  value  for  the  free  molecule  shear  stress  is  given 
Incorrectly  hy  Equation  (8)  -  it  is  approximately  15  percent  too  large. 

With  this  correction  the  shear  stress  is  seen  to  he  in  good  agreement  with 
the  Reference  2  value  at  large  d/X. 

In  sumsary,  the  simple  discrete  velocity  model  gives  a  reasonably  accu¬ 
rate  description  of  low  Mach  msnber  Couette  flow  over  the  complete  range  of 
d/X.  It  should  he  remembered,  however,  that  the  calculations  of  the  col¬ 
lision  cross  sections  in  the  Appendix  can  only  he  defended  as  being  reasonable. 
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A  further  study  of  this  flow  with  more  molecular  velocity  cells  is  required 
before  we  can  he  sure  that  the  accuracy  of  the  eight  cell  solution  is  not 
simply  fortuitous. 


1.0 


d/X 


Fl^ur*  4. 

Velocity  slip  and  ahear  rtrta*  m  *  function  of  BoudMC  mater 
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AFPKHDn 


To  find  the  rate  of  change  of  due  to  colUalons  we  oust  first  de¬ 
termine  which  of  its  possible  collisions  cause  a  net  loss  or  gain  front  cell 
1  as  distinguished  frost  those  in  which  the  collision  partners  simply  exchange 
cells.  Refer  to  Table  1  in  which  the  sign  of  the  velocity  components  are 
listed.  Molecules  nay  Jump  only  between  pairs  of  cells  In  which  there  is 


>cell 

1 

2 

3 

4 

5 

6 

7 

8 

u 

- 

+ 

- 

+ 

- 

+ 

- 

+ 

V 

+ 

+ 

m 

+ 

+ 

m 

m 

w 

+ 

+ 

+ 

+ 

m 

- 

- 

m 

Table  1.  Sign  of  velocity  components 

the  same  combination  of  signs  for  u,  v,  and  w.  Thus  the  collision  1-4  -*-2-3 
is  possible  but  the  collision  1-2  has  no  effect.. 

The  collisions  in  which  is  depleted  are: 

1-4 -.2-3  1-6 -.2- 5  1-7  -3-5 
1-8 -.2- 7  1-8 -.3-6  i-8— 4-5 

and  the  rate  of  loss,  is 

\  -  vr‘  8/  (Hx  *4  +  *!  »6  +  »!  V  +  v\  C  (H1  V 
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where  the  superscript  a  refers  to  the  relative  velocity  and  effective  cross 
section  of  collision  partners  diagonally  opposed  on  the  faces  of  the  velocity 
Cube  In  Pig .  2  and  b  to  those  on  the  cube  diagonals. 

Fran  synwtry,  molecules  from  a  collision  1-8,  for  Instance,  are  equally 
distributed  to  2-7,  3-6,  and  4-5,  and  therefore  the  rate  of  addition  to  cell 

1  is 

°!  ■  V  S>2  "3  *  "2  V  *3  V 
♦  j  Tpk  s.1  (h2  s,  +  »3  »6  ♦  n,  »5) 

The  relative  velocities  are  given  by 

Yr‘  ■  2/l  S  ^  ■  2  «• 

The  effective  collision  cross  sections  are  determined  for  hard  elastic* 
spheres  as  follows.  Consider  first  an  encounter  of  the  kind  1-4,  which 
takes  place  In  the  v  -  constant  plane,  as  sketched  In  Pig.  lb.  The  question 
is  what  fraction  of  those  molecules  that  collide  should  be  assigned  to  cells 

2  and  3  (with  the  others  being  returned  to  1  and  4).  Since  the  scattering 
has  circular  symmetry  the  most  reasonable  assumption  is  that  half  of  the 
colliding  molecules  are  deflected,  l.e.,  that 

3.*  -  i  3. 

where  S  is  the  mutual  collision  cross  section.  The  corresponding  assumption 
for  collisions  such  as  1-6,  In  which  there  Is  scattering  in  three  dimensions, 
leads  to  the  result 


Seb  -  1 4/T  S. 

Finally,  waking  use  of  the  syametry  about  w  ■  0,  i.e.,  of  the  condition 
N1  ■  Hg  -  Mg,  etc.,  and  dividing  and  1*^  by  N,  ve  get 

<ST>e  '  5  <°i  '  V  ’  (1  +  “  “<”2  V  "l  "t 

-  2  9  (n2  n3  •  nx  n^) 

where 

e-~.(i+|.|cSN  M 

4*  t  fil  W 

The  solutions  of  Eqs.  (5a)  and  (5b)  or  Eqs.  (l8a)  and  (l8c)  of  Ref.  1 
and  hence  quantities  such  as  viscosity,  velocity  slip,  etc.,  are  given 
directly  In  terns  of  6  but  the  fluid  mechanical  results  are  usually  stated 
in  terns  of  the  mean  free  path,  1.  To  introduce  this  parameter  ve  have  to 
assume  that  the  gas  is  close  to  equlllbrluB  and  that  the  velocity  gradients 
are  snail,  l.e.,  that  in  our  case  the  M^'s  are  nearly  equal. 

The  collision  rate  of  all  the  molecules  of  cell  1,  for  Instance,  is 
given  by 

9  -  e  8  1 2/fi (^  Hg  +  1^  R3  +  ^  M?)  +  2  fijl  (^  ^  Mg  +  ^ 

♦  2  >1  »$] 

-  (6/^3  +  ^2/3  +  2)  c  8  l^2  -  10.37  c  8  Mx2 


-16. 


I 


Then  the  collision  frequency  per  molecule,  o,  ie 


a 


10.37  cBIIj*  8  c 


end 

x  .£  „ _ 1_ 

*  o  1,3  81' 

Nov  recall  that 

8  -  ~=  (1  +  2//1)  c  S  N  -  .88  c  8 

r 


.68  o  -  .68  f  . 


and  therefore 


1.3 


a  - 
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H  -  1.3  c  8  H, 


1 

SYMBOLS 

c 

aeon  thermal  speed 

d 

distance  between  vails  In  Couette  flow 

f 

distribution  function 

a 

gain  in  f  due  to  collisions 

H 

Boltzmann  H  function 

L 

loss  in  f  due  to  collisions 

m 

molecular  weight 

R 

number  of  molecules  per  unit  voluse 

"i 

lumber  of  molecules  per  unit  rolime  with  velocity  ^ 

nl 

fraction  of  molecules  per  unit  voluse  with  velocity  ^ 

pyx 

shear  stress 

*pyx^f.B. 

<1 

free  molecule  shear  stress 

c//T,  component  of  c  along  coordinate  axes 

8 

mutual  collision  cross  section 

8e 

effective  collision  cross  section 

t 

time 

U 

x  component  of  fluid  velocity 

difference  in  vail  velocities 

AU 

velocity  slip 

u 

x  component  of  molecular  velocity 
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SYMBOLS  (Cont'd.) 


<P 

X 

V> 

P 

a 


y  component  of  fluid  velocity 
y  component  of  molecular  velocity 
molecular  velocity 
relative  molecular  velocity 
b  component  of  molecular  velocity 
coordinate  axes 
y/d 

d  e/2  q 

effective  colllalon  rate, 

collision  rate  of  occupants  of  one  velocity  cell,  per  unit 
volume 

mean  free  path 
viscosity 
density,  m  H 

collision  frequency  per  molecule  ,  c/X 


■==■  (l  +  3/^*3)  c  S  H  »  .68  c/X 

•J6 


Superscripts 

a  refers  to  collision  partners  on  sane  face  of  velocity  cube, 

Figure  2 

b  refers  to  collision  partners  on  velocity  cube  diagonals. 

Figure  2 
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